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This  Is  the  first  part  of  an  Investigation  dealing  with  the  global  load- 
displacement  behavior  of  a  stationery  center  cracked  panel  in  which  materiel 
damage  Is  accumulated  for  each  Increment  of  loading.  The  finite  element  pro¬ 
cedure  was  used  where  the  local  stiffness  In  regions  ahead  of  the  crack  under¬ 
go  changes  In  accordance  with  a  predetermined  failure  criterion.  The  result¬ 
ing  non-linear  behavior  of  the  load-displacement  curves  are  obtained  by  two 
different  damage  models.  Growth  of  the  crack  will  be  considered  In  the  second 
rt  of  this  work.  -\ 
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The  proposed  damage  concept  models  the  development  of  microcracks  and/or 
voids  ahead  of  a  macrocrack  by  altering  the  local  stiffness.  As  the  material 
elements  are  damaged,  the  linearly  elastic  behavior  changes  to  pseudo- lin¬ 
early  elastic  behavior,  under  the  assumption  that  any  energy  dissipated  due 
to  microcrack  generation  Is  no  longer  recoverable  upon  unloading.  No  other 
dissipative  mechanism  Is  assumed,  j 


-A  Single  Damage  Level  model  was  based  on  a  "damaged/undamaged"  evaluation 
of  each, finite  element  In  the  panel  specimen.  The  damage  threshold  was  based 
on  the  elastic  strain  energy  density  At  the  yield  stress  level  of  an  uniaxial 
tensile  test  specimen.  The  effect  of  damage  on  the  bulk  material  properties 
In  the  damage  zone  were  taken  to  be  a  reduced  secant  elastic  modulus  consis¬ 
tent  with  the  uniaxial  tensile  test  specimen  and  with  the  effective  elastic 
properties  of  the  mlcrocracked  portion  of  the  medium. .  The  damage  status  of 
the  element  was  the  Independent  variable  of  the  modeu  subject  to  the  fail¬ 
ure  sequence  of  all  elements  determined  by  the  strain  energy  density  ranking 
of  each  element.  The  applied  load  necessary  former given  element  to  reach 
the  strain  energy  density  value  at  yield  stress  was  the  dependent  quantity 
of  the  model. _ _ — " 

£  A  Multiple  Damage  Level  model  Incorporated  twenty-four  discrete  values 
of  the  secant  elastic  modulus  of  a  specified  uniaxial  teiplle  specimen  true 
stress-true  strain  curve.  The  strain  energy  density  level  at  yield  stress 
was  taken  to  be  the  damage  threshold  level. w^The  value  of  the  secanlr modulus 
representing  the  highest  level  of  material  duagejsas  chosen  to  be  consistent 
with  the  ultimate  stress  and  ultimate  strali  \valueso?  the  tensile  specimen 
curve,  and  the  effective  elastic  properties  of  the  mlcrocracked  portion  of 
the  medium.  The  applied  load  was, the  Independent  variable  which  determined 
the  level  of  damage  In  each  element  of  the/model  for  a  given  load.  Note 
that  this  Is  In  contrast  to  the  Single  Dange  Level  model,  where  the  reverse 
was  true.  •  \ 
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ABSTRACT 


This  Is  the  first  part  of  an  Investigation  dealing  with  the  global  load- 
displacement  behavior  of  a  stationary  center  cracked  panel  In  which  material 
damage  Is  accumulated  for  each  Increment  of  loading.  The  finite  element  proce¬ 
dure  was  used  where  the  local  stiffness  In  regions  ahead  of  the  crack  undergo 
changes  In  accordance  with  a  predetermined  failure  criterion.  The  resulting 
non-linear  behavior  of  the  load-displacement  curves  are  obtained  by  two  differ¬ 
ent  damage  models.  Growth  of  the  crack  will  be  considered  In  the  second  part 
of  this  work. 

The  proposed  damage  concept  models  the  development  of  microcracks  and/or 
voids  ahead  of  a  macrocrack  by  altering  the  local  stiffness.  As  the  material 
elements  are  damaged,  the  linearly  elastic  behavior  changes  to  pseudo-11 nearly 
elastic  behavior,  under  the  assumption  that  any  energy  dissipated  due  to  micro¬ 
crack  generation  Is  no  longer  recoverable  upon  unloading.  No  other  dissipative 
mechanism  Is  assumed. 

A  Single  Damage  Level  model  was  based  on  a  "damaged/undamaged"  evaluation 
of  each  finite  element  In  the  panel  specimen.  The  damage  threshold  was  based 
on  the  elastic  strain  energy  density  at  the  yield  stress  level  of  an  uniaxial 
tensile  test  specimen.  The  effect  of  damage  on  the  bulk  material  properties  In 
the  damage  zone  were  taken  to  be  a  reduced  secant  elastic  modulus  consistent  with 
the  uniaxial  tensile  test  specimen  and  with  the  effective  elastic  properties  of 
the  mlcrocracked  portion  of  the  medium.  The  damage  status  of  the  element  was 
the  Independent  variable  of  the  model,  subject  to  the  failure  sequence  of  all 
elements  determined  by  the  strain  energy  density  ranking  of  each  element.  The 
applied  load  necessary  for  a  given  element  to  reach  the  strain  energy  density 
value  at  yield  stress  was  the  dependent  quantity  of  the  model. 
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A  Multiple  Damage  Level  model  Incorporated  twenty-four  discrete  values  of 
the  secant  elastic  modulus  of  a  specified  uniaxial  tensile  specimen  true  stress- 
true  strain  curve.  The  strain  energy  density  level  at  yield  stress  was  taken  to 
be  the  damage  threshold  level.  The  value  of  the  secant  modulus  representing  the 
highest  level  of  material  damage  was  chosen  to  be  consistent  with  the  ultimate 
stress  and  ultimate  strain  values  of  the  tensile  specimen  curve,  and  the  effec¬ 
tive  elastic  properties  of  the  mlcrocracked  portion  of  the  medium.  The  applied 
load  was  the  Independent  variable  which  determined  the  level  of  damage  In  each 
element  of  the  model  for  a  given  load.  Note  that  this  Is  In  contrast  to  the 
Single  Damage  Level  model,  where  the  reverse  was  true. 
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INTRODUCTION 


There  Is  Increasing  concern  In  the  field  of  continuum  mechanics  that  the 
Interaction  between  changes  In  material  properties  and  material  behavior,  due 
to  damage  In  the  form  of  microcracks  or  other  defects.  Is  not  being  distinguished 
and  assessed  quantitatively.  In  recent  years,  attention  has  been  focused  on 
studying  the  accumulation  of  physical  defects  in  the  solid  as  It  Is  being  loaded 
and  observing  the  changes  In  the  state  of  stresses  and  strains.  This  differs 
fundamental 1y  from  the  classical  formulation  of  constitutive  equations  that  do 
not  address  the  physical  damage  of  the  material. 

In  order  to  develop  more  realistic  constitutive  relations  that  reflect  the 
occurrence  of  Irreversibility  caused  by  physical  damage  of  the  loaded  material, 
the  behavior  of  a  center  cracked  panel  (Figure  1)  under  unidirectional  tensile 
loading  perpendicular  to  the  crack  has  been  studied  In  detail.  Such  a  specimen 
features  a  macro-defect  In  the  form  of  a  symmetrically  located  crack,  whose 
length  Is  of  the  same  order  of  magnitude  as  the  specimen  size.  Under  sufficient¬ 
ly  high  loads,  the  macrocrack  would  be  expected  to  Increase  In  length.  For  some 
value  of  the  load,  the  process  will  become  unstable  and  result  In  failure  of  the 
panel . 

A  specimen  made  of  a  "ductile"  material  would  be  expected  to  develop  damage 
at  a  lower  scale  level  In  the  form  of  microcracks  and/or  microvoids.  These  de¬ 
fects  may  be  Initially  present  In  the  material,  or  may  coalesce  from  the  move¬ 
ment  of  dislocations  and  other  defects  at  an  even  smaller  scale,  l.e.,  the  atomic 
scale.  In  the  presence  of  a  macrocrack,  such  regions  are  known  to  exist  ahead 
of  the  crack  tip  on  both  sides  of  a  two-ended  crack.  The  latter  regions,  common¬ 
ly  Identified  as  "plastic  enclaves"  in  the  classical  theory,  of  plasticity,  have 
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Figure  1  -  Center  cracked  panel  specimen 


been  the  subject  of  much  attention  since  the  dimensions  of  such  regions  can  be 

significant  when  compared  to  the  crack  length.  The  microcracks  themselves  will 

produce  microplastic  regions  ahead  of  both  crack  tips,  also.  Their  effect  on 

★ 

the  propagation  of  the  macrocrack  Is  not  yet  fully  understood  . 

In  order  to  develop  a  model  of  the  processes  described  above.  It  Is  assumed 
that  the  strain  energy  density,  dkl/dV,  per  unit  volume  of  material.  Is  the  quan¬ 
tity  that  reflects  material  damage  at  the  micro  level,  as  well  as  the  length  and 
direction  of  macrocrack  propagation  at  the  macro  level**.  The  term  "material  dam¬ 
age”  will  refer  to  processes  which  utilize  strain  energy  at  the  continuum  scale  to 
produce  Irreversible  changes  In  the  material.  These  changes  will  In  turn  affect 
the  local  material  behavior  as  a  result  of  the  formation  of  microcracks,  micro¬ 
voids  or  coalescence  of  these  defects. 


Delameter,  Herrmann  and  Barnett  have  considered  the  two  dimensional  problem 
of  a  doubly  periodic  array  of  cracks  with  equal  length  2a.  The  cracks  are  paral 
lei  to  the  x-axIs  and  perpendicular  to  the  y-axls.  Dislocations  have  been  used 
to  describe  the  discontinuity  In  the  displacement  field  which  exists  across  the 
cracks.  The  resulting  singular  Integral  equations  of  the  first  klnc  have  been 


The  terminology  "small  scale  yielding"  refers  to  relative  size  of  the  plastic 
enclave  as  compared  with  crack  length  and  has  caused  much  confusion  when  using 
elastic-plastic  solutions  of  a  crack  In  an  Infinite  medium,  l.e.,  containing  no 
specimen  dimensions.  It  has  been  shown  In  [9]  and  discussed  In  [13]  that  the 
plastic  enclaves  can  be  many  times  greater  than  the  crack  length  and  the  state 
of  affairs  near  the  crack  Is  still  essentially  elastic  because  the  stored  elas¬ 
tic  energy  In  an  Infinite  medium  ahead  of  the  crack  Is  always  much  larger  than 
the  stored  elastic-plastic  energy  In  the  enclaves  that  are  finite  In  dimension. 

An  additional  length  parameter  other  than  the  crack  length  such  as  specimen  dimen¬ 
sion  or  plate  thickness  Is  required  for  analyzing  the  Influence  of  plasticity  on 
the  global  behavior  of  a  specimen  or  structural  component. 

*• 

It  Is  considered  to  be  pertinent  that  the  failure  criterion  at  both  the  micro 
and  macro  level  be  consistent  In  order  to  minimize  the  Introduction  of  arbitrari¬ 
ness  Into  the  model. 
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numerically  solved.  This  formulation  accounts  for  the  Interaction  of  the  periodic 
cracks.  The  effective  elastic  constants  of  the  medium  with  elastic  modulus  E, 
shear  modulus  G,  and  Poisson's  ratio  v  are  calculated  on  the  basis  of  the  strain 
energy  density  per  unit  thickness.  In  plane  strain,  the  relation  for  strain  en¬ 
ergy  density  Is 


dW 

W 


x  + 


(1) 


for  a  medium  In  the  absence  of  cracks,  and  Is  equated  to  the  relation  for  a 
cracked  medium. 


v12°xay  . 
F  + 


(2) 


Here,  AE/bd  Is  the  energy  density  change  due  to  one  crack  In  the  material,  with 
a  spacing  between  adjacent  crack  centers  of  b  In  the  y-dlrectlon  and  d  In  the  x- 
dlrectlon.  For  a  crack  subject  to  Mode  I  loading, 

AEj  »  xa2a2(l-v2)BI/E  (3a) 


while  for  Mode  II  loading. 


AEn  •  *e2a2(l-v2)BII/E 


(3b) 


where  B1  and  B11  are  coefficients  which  are  functions  of  the  distances  b  and  d 
between  cracks.  For  mixed  mode  loading,  the  total  energy  Is  given  as 

AE  -  AEj  +  AEn 
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The  bulk  properties  of  the  cracked  medium  are  orthotropic  and  calculated  to  be: 


E 


x 


E 


fr.Ill _ '  _ 

E  v  l+(2ira2BVbd) 

6xv  _ _ 1 

*  H-te^B^G/bdE) 

The  asymptotic  values  of  Ey  and  are  found  to  be 
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(2a/b)-~ 


VE 


and 


(4a) 

(4b) 

(4c) 

(4d) 


(5a) 
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(2a/b)-~ 
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(5b) 


The  first  expression  Implies  that  a  value  of  E^/E  greater  than  0.5  will  be  the 
asymptotic  value  of  E^/E  unless  the  uncracked  ligament  d-2a  between  any  two 
colllnear  cracks  becomes  as  small  as  the  crack  periodicity,  b.  In  the  y-dlrec- 
tlon.  This  Is  one  Indication  that  high  crack  densities  may  still  result  In  a 
significant  effective  elastic  modulus  of  the  bulk  material.  The  shear  modulus 
Is  reduced  significantly  by  this  formulation,  to  a  limiting  value  of  zero. 


Gottesman,  Hashln  and  Brull  have  Investigated  the  upper  and  lower  bounds  of 
the  effective  properties  of  a  cracked  medium  using  the  principals  of  minimum  po- 


tential  and  minimum  complementary  energy,  respectively.  For  two  dimensional 
plane  strain,  the  change  In  energy  of  the  total  body  under  uniform  loading  In 
the  presence  of  a  statistically  homogeneous  crack  distribution  will  be  for  m 
total  cracks. 


al  V19°v°u  TSu 

9  m  "  *  k  ‘  <2F  *  &  '  "HT1  +  T*>v  *  i  % 


(6) 


where  V  Is  the  volume  of  the  body  and  AWm  Is  the  energy  change  associated  with 
each  crack.  For  a  crack. 


“V  ■  T  "(J  /  C“iVs  <7) 

Tn 
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where  [u^]  Is  the  displacement  discontinuity  across  the  crack  face,  Sm  Is  the 
nth  crack  surface  and  nj  Is  the  unit  normal  to  the  crack  surface.  Alternately, 
In  terms  of  the  stress  intensity  factor,  the  following  expression  Is  obtained: 


where  1$  a  crack  length  parameter  and  KI’  KII  and  Kjjf  are  the  Mode  I,  II  and 
III  stress  Intensity  factors,  respectively. 

If  the  uncracked  body  Is  Isotropic,  and  If  the  crack  distribution  Is  sta¬ 
tistically  Isotropic  also,  then  the  effective  bulk  modulus  IT  and  shear  modulus 
5  can  be  expressed  In  the  form 


and 


l 

G 


1  +  1 

G  +  IG 


(9b) 


where  the  quantities  1/aK  and  1/AS  are  functions  of  the  displacement  disconti¬ 
nuities  across  the  crack  surface  and  the  globally  applied  stress. 

For  the  case  of  low  crack  density,  the  upper  bound  calculation  is  performed 
using  the  principal  of  minimum  potential  energy.  This  solution  does  not  con¬ 
sider  crack  Interaction.  A  crack  density  parameter  a,  of  the  form 


a 


(10a) 


for  two  dimensional  problems,  and 


«4e  a  b2  (10b) 

V  m  mm 

for  three  dimensional  problems  are  Introduced,  where  am  and  bm  are  crack  dimen¬ 
sion  parameters.  The  results  for  randomly  oriented  circular  cracks  are 

Eupper  *  E/Cl  +  1M1-v2)(10-3v)/45(2-v)]  (11a) 

6upper  *  W  *  32a(1-v)(5-v)/45(2-v)]  (lib) 

where,  for  circular  cracks  of  radius  am 

TO 
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The  lower  bounds  are  formulated  In  terms  of  the  material  properties  of  the 
uncracked  medium,  the  crack  density  parameter  a,  and  crack  geometry  parameters. 
This  last  condition  Is  quite  restrictive  In  calculating  lower  bounds  since  It 
requires  Information  on  each  crack  In  the  general  case  of  random  geometry  cracks, 
l.e.. 


Elower  =  E  /  f^a*am‘bm’cm^dam 

where  am,  bm  and  cm  are  crack  geometry  parameters. 

Despite  the  difficulties  Inherent  In  the  lower  bound  formulation,  the  upper 
and  lower  bounds  for  selected  crack  geometries  are  found  to  be  within  a  few  per¬ 
cent  of  each  other  In  most  cases  considered.  The  results  found  by  Delameter, 
Herrmann  and  Barnett  lie  within  the  predicted  upper  and  lower  bounds. 

Budlansky,  O'Connell  and  Hoenlg  have  applied  the  self-consistent  procedure 
to  the  problem  of  dry,  as  well  as  moisture  saturated  cracks  In  an  elastic  medium. 
By  equating  the  potential  energy  j  In  a  representative  volume  element  of  a  solid 
with  distributed  cracks  to  the  potential  energy  +  In  the  original  material  plus 
the  change  a*  due  to  the  cracks,  an  expression  of  the  form  Is  obtained: 

03) 

2K  I  m 

where  the  quantities  with  the  over  bar  refer  to  the  cracked  medium.  In  equation 
(13),  p  Is  a  hydrostatic  state  of  stress,  V  Is  the  volume  of  the  body,  K  Is  the 
bulk  modulus,  E  Is  the  elastic  modulus,  v  Is  Poisson's  ratio,  and  a  Is  a  crack 
length  parameter.  The  function  f  Is  a  non-dimensional  shape  factor  dependent  on 
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crack  geometry  . 


For  the  case  of  uniaxial  tension  aQ  applied  to  the  body. 


-  7T~  Cz1  £  a3f(v)  +  ^  z  a3g(v,6)] 

2^  “  £  IB  f  ID 


Another  non-dimensional  function  g  Is  Introduced  In  conjunction  with  the  resolved 
local  shear  stress  x  at  the  crack  edge.  The  angle  a  Is  dependent  on  the  orienta¬ 
tion  of  the  cracks.  The  function  f  Is  associated  with  the  resolved  local  tensile 
stress  o  perpendicular  to  the  crack  edge. 

The  effective  elastic  modulus  ?,  effective  bulk  modulus  7,  and  effective 
Poisson's  ratio  v  are  assumed  to  follow  the  conventional  relation 


r  -  3l((l-2v) 


Thus,  three  equations  and  three  unknowns  exist  to  determine  the  material  con¬ 
stants  of  the  cracked  medium.  The  effective  shear  modulus  £  can  be  calculated 
from  these  three.  Assuming  that  some  knowledge  of  crack  size  and  geometry  Is 
known.  In  the  statistical  sense,  the  summation  terms  In  the  potential  energy 
expressions  can  be  replaced  by  correlated  or  uncorrelated  averages  of  a,  f(v), 
and  g(v,e). 

In  order  to  evaluate  the  functions  f  and  g,  the  average  crack  energy  Is 
calculated  since  this  Is  the  negative  of  the  potential  energy  change.  The  ex- 


Note  that  the  critical  assumption  made  Is  that  the  crack  behaves  as  If  It  lies 
In  the  medium  Influenced  by  the  other  cracks.  This  medium  exhibits  the  effective 
material  properties  of  the  cracked  medium.  Thus,  crack  Interaction  Is  addressed, 
but  without  considering  each  crack  Individually. 


V 


presslon  discussed  above,  equation  (8),  which  utilizes  the  three  stress  Intensity 
factors.  Is  used  for  these  calculations.  Once  this  Is  accomplished  for  the 
crack  geometries  under  consideration,  the  material  properties  can  be  calculated 
using  a  crack  density  parameter 


e 


2N 

it 


<£> 


(16) 


where  N  Is  the  total  number  of  cracks,  A  Is  the  crack  area,  P  the  crack  perim¬ 
eter,  and  the  brackets  denote  ensemble  averaging.  This  Is  similar  to  the  one 
cited  by  Gottesman,  Hashln  and  Brull,  as  can  be  seen  for  the  case  of  all  circular 
cracks  In  the  medium,  for  which 


e  *  N  <a3> 
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The  self-consistent  method.  Interestingly  enough,  predicts  a  value  of  c  9/16 
which  results  In  v  +  0,  which  Is  shown  to  Imply  F/E  and  tyK  ♦  0.  Thus,  while 
not  a  rigorous  criteria  for  material  failure,  this  value  of  e  could  be  Inter¬ 
preted  as  corresponding  to  a  condition  of  Instability  In  the  material  based  on 
crack  size  and  population. 

For  all  values  of  0^v<p.5,  the  material  properties  are  found  to  be  preserved 
to  at  least  40%  of  their  original  value  If  e  does  not  exceed  0.2.  This  Is  for 
the  case  of  dry  (moisture  unsaturated)  randomly  distributed  circular  cracks. 

The  behavior  of  a  medium  with  a  distribution  of  cracks  can  be  formulated  In 
terms  of  what  Elmer  has  labelled  as  pseudo-linear  elasticity.  Recall  that  a 
defect  free  medium  Is  linearly  elastic  If  It  Is  poi/JUvtty  homogtnttnu,  l.e.. 
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and,  moreover,  additive. ,  l.e. , 

V*!]’  ♦  4j2)l  ■  *  •ijt«g)3  ('“>) 

where  the  superscripts  1  and  2  denote  arbitrary  admissible  strain  fields  for  a 
given  body  and  boundary  conditions.  If  the  medium  contains  defects  such  as  mi¬ 
crocracks,  the  additive  condition  may  be  assumed  to  be  violated.  This  violation 

Is  consistent  with  the  Introduction  of  an  strain  field  which  opens  the  cracks 

(21 

In  a  portion  of  the  body,  followed  by  the  superposition  of  a  strain  field  cjj' 
which  would  close  some  or  all  of  the  cracks  which  were  open.  If  the  condition 
of  positive  homogeneity  Is  retained  (Implying  that  open  cracks  would  remain  open 
and  closed  cracks  would  remain  closed  upon  amplification  or  attenuation  of 
by  a  positive  factor  k)*the  material  behaves  In  a  pseudo-linear  fashion.  The 
crack  closure  process  Is  considered  reversible  on  the  assumption  that  no  tangen¬ 
tial  forces  develop  during  closure.  The  material  behavior  Is  thus  reversible 
upon  unloading  from  both  tensile  and  compressive  stress  fields.  Only  microcrack 
propagation  may  be  assumed  to  dissipate  elastic  energy  since  the  material  Itself 
Is  by  definition  linear. 

In  this  Investigation,  the  results  for  a  center  cracked  panel  specimen  In  the 
absence  of  macrocrack  propagation  are  presented.  Macrocrack  growth  Is  being  con¬ 
sidered  and  will  be  reported  separately.  The  Strain  Energy  Density  Theory  of 
fracture,  as  proposed  by  Slh  [12,13]  will  be  used  to  predict  macrocrack  growth 
In  conjunction  with  the  microdamage  models  to  be  discussed  below.  This  approach 
to  the  problem  separates  the  effects  of  micro  and  macrodamage.  By  comparing  these 
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results  with  the  combined  effects  of  damage  at  both  levels,  understanding  Is 
gained  as  to  which  level  of  damage  predominates  the  specimen  global  failure 
process  for  different  structural  materials. 

The  primary  concerns  In  this  Investigation  have  been  (1)  the  Influence  of 
damage  on  effective  material  properties,  and  the  relation  of  this  to  (2)  the 
global  specimen  behavior.  Not  directly  addressed  at  this  point,  but  of  parallel 
concern,  are  the  means  by  which  the  effective  properties  of  a  material  In  both 
the  virgin  and  damaged  states  are  obtained  experimentally.  These  and  subsequent 
analyses  should  provide  greater  understanding  of  the  relation  of  tensile  test 
specimen  data  to  global  structural  behavior. 


DAMAGE  MODELS  AND  FINITE  ELEMENT  FORMULATION 


Generally  speaking,  material  nonlinearity*  exhibited  on  the  stress-strain 
curve  for  metals  Implies  Irreversibility  due  to  material  damage.  This  being  In¬ 
herently  a  path  dependent  process  Is  sensitive  to  the  rate  of  loading  as  each 
Increment  of  load  Increase  will  affect  the  stress  and  strain  state  which  In 
turn  governs  the  next  Increment  of  material  damage.  Hence,  nonlinear  stress 
and  failure  analysis  cannot  be  analyzed  Independently  as  In  the  linear  case  and 
must  be  considered  In  sequence  within  each  Increment  of  loading.  Failure  to  do 
so  results  In  unphysical  and  unrealistic  prediction. 

All  stress  analysis  used  to  generate  the  load-displacement  curves  were  car¬ 
ried  out  on  the  center  cracked  panel  specimen  for  a  material  with  characteristics 
of  a  low  yield  stress  (<iy  •  35,000  psl)  structural  steel.  The  load-displacement 
curves  were  calculated  Incrementally  In  a  series  of  "stress  analysis-material 
damage  assessment"  Increments  or  cycles. 

The  Axlsymnetrlc/Planar  Elastic  Structures  (APES)  two-dimensional  fracture 
mechanics  and  stress  analysis  finite  element  program  [8]  performed  the  stress 
analyses  portion  of  each  Increment.  The  formulation  of  APES  utilizes  quad-12 
elements  which  allow  for  cubic  displacement  fields  and  quadratic  stress  and 
strain  fields  within  each  element.  The  r^2  displacement  field  In  the  Immediate 
vicinity  of  the  crack  tip  Is  embedded  In  the  solution  through  the  use  of  1/9  to 
4/9  nodal  spacing  on  the  element  sides  adjacent  to  the  crack  tip.  The  tensile 
loading  of  the  panel  subjects  the  crack  to  a  symmetric  Mode  I  situation,  thus 
requiring  only  one  quarter  of  the  physical  problem  (Figure  2)  to  be  modelled  by 

Global  and  local  nonlinearity  should  be  distinguished.  The  present  analysis 
assumes  linearity  for  the  undamaged  material  while  nonlinear  response  Is  devel¬ 
oped  globally  for  each  Increment  of  loading. 
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Figure  2  -  Finite  element  grid 


440  nodes  which  define  77  elements.  The  CDC  6400  computer  used  In  this  study 
typically  required  225  system  seconds  and  132k  of  central  memory  to  complete 
the  stress  analysis  for  a  given  specimen  load  and  damage  distribution. 

The  damage  processes  that  occur  In  the  material  manifest  themselves  by 
changes  In  the  material  properties  at  the  continuum  level.  The  finite  element 
formulation  allows  for  this  by  associating  with  each  element  one  of  up  to  25  dif¬ 
ferent  material  property  pairs  (l.e.,  elastic  modulus  E,  and  Poisson's  ratio  v). 
Thus,  a  set  of  properties  (E^.v-j)  corresponding  to  the  Initial  elastic  behavior 
of  the  material,  as  defined  from  an  uniaxial  tensile  test  specimen,  and  up  to  24 
discrete  material  property  pairs  (E^.v^  for  1  =>  2, 3,4,..., 24,  different  Instances 
In  the  damped  material's  developing  history. 

Two  distinct  types  of  material  damage  models  were  used  to  generate  load-dis¬ 
placement  curves  for  the  center  cracked  specimen.  The  exact  type  of  damage  which 
accumulates  In  the  material,  l.e.,  the  microcrack  sizes  and  shapes,  has  not  been 
specified  In  detail.  Rather,  reasonable  lower  bounds  on  the  effective  properties 
of  the  damaged  material  have  been  used  to  generate  the  bulk  behavior  of  cracked 
elastic  media.  The  type  of  damage  producing  this  bulk  behavior  obviously  Is  not 
unique.  In  the  sense  that  a  wide  variety  of  crack  distributions,  geometries, 
sizes,  and  orientation  can  result  In  the  same  behavior  of  the  cracked  media.  The 
rate  of  loading  can  also  alter  the  bulk  behavior. 

The  first  type  of  model  provided  one  of  two  conditions  for  each  element  in 
the  center  cracked  panel  specimen,  l.e.,  either  "undamaged"  or  "damaged".  This 
simple  model  is  Interesting  In  two  respects.  First,  the  single  level  of  damage 
reduces  the  complexity  of  the  developing  damage  zone.  This  allows  the  effect 
of  damage  In  a  single  region  (or  element  In  this  case)  to  be  directly  evaluated 
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without  additional  changes  In  material  properties  In  other  regions  of  the  panel 
specimen.  Second,  the  damage  process  Is  considered*  as  a  sequence  of  "events", 
l.e.,  changes  In  element  material  properties  from  the  undamaged  to  the  damaged 
states.  Thus,  the  damage  process  Is  discretized  to  one  element  at  a  time. 


The  damage  In  the  panel  specimen  was  assessed  In  the  method  outlined  below. 
After  each  stress  analysis  of  the  center  cracked  panel,  the  undamaged  elements 
were  ordered  In  terms  of  decreasing  average  strain  energy  density.  Ideally, 
each  element  experiencing  damage  would  specify  an  applied  load  Increment,  thus 
producing  the  unique  element  damage  sequence  and  load-displacement  curve  for  each 
panel  specimen.  Due  to  constraints  on  computing,  however,  multiple  element  dam¬ 
age  was  used  to  generate  the  load  Increment  subject  to  the  upper  and  lower  bound¬ 
ing  criteria  for  the  load  Increments  discussed  In  Appendix  A.  By  Incrementing 
the  damage  of  the  elements  beyond  the  Initial  Increment  In  groups  typically  no 
longer  than  five  elements  per  Increment,  load  displacement  curves  consistent  with 
the  damage  model  and  bounding  procedures  were  obtained.  The  change  aS/Vj  In  panel 
stiffness  S  per  the  1th  damaged  element  volume  was  small  when  compared  to  the 
overall  panel  stiffness  value.  In  the  event  that  the  damage  region  of  the  pre¬ 
vious  increment  enhanced  the  dW/dV  value  of  an  undamaged  element  to  a  value  Indi¬ 
cating  failure  at  the  prior  load,  the  failure  of  such  an  element  could  be  Incor¬ 
porated  Into  the  next  damage  Increment  with  a  negligible  effect  on  the  overall 
characteristics  of  the  load-displacement  curve. 

Two  analyses  of  the  center  cracked  panel  were  performed  using  the  Single  Dam¬ 
age  Level  model.  In  both  cases,  the  strain  energy  density  at  the  yield  stress, 
(dW/dV)y$  •  o*/2E,  of  a  uniaxial  tensile  specimen  was  used  as  the  criterion  by 


To  the  ordinary  mind,  the  material  damage  Is  a  discrete  process  for  the  simple 
reason  that  continuity  cannot  be  conceived  physically  but  only  In  the  limit  as  a 
mathematical  process. . 


.1  - 
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which  the  damage  process  would  occur  In  the  material.  This  Is  consistent  with 
the  assumption  that  the  nonlinear  behavior  of  the  tensile  specimen  coincides  with 
the  onset  of  elastic  energy  dissipation  to  damage  processes.  The  damage  status 
of  the  1th  element  Is  determined  according  to  the  following 


({jjjr)  1  (37)  Undamaged  Material 


ys 


^  >  (W}ys  Damaged 


Material 


09a) 

(19b) 


In  the  first  application  of  the  model,  the  damaged  material  elastic  modulus 
was  assumed  to  be  a  factor  of  0.8  times  the  undamaged  modulus  obtained  from  the 
uniaxial  tensile  specimen  (Figure  3).  For  the  second  application  of  the  model, 
the  damaged  material  elastic  modulus  was  taken  to  be  a  factor  of  0.6  times  the 
undamaged  modulus  (Figure  4).  Since  crack  propagation  Is  not  addressed  In  these 
analyses,  the  exact  form  of  the  material's  true  stress-true  strain  curve  (Figures 
3,4)  Is  rather  arbitrary  beyond  the  yield  point.  This  is  due  to  the  Independence 
of  the  microdamage  model  on  the  critical  strain  energy  density,  (dW/dV)cr,  which 
Is  attained  at  fracture  of  the  tensile  specimen. 
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The  second  damage  model  used  to  generate  a  load-displacement  curve  for  the 
center  cracked  panel  specimen  was  based  on  24  discrete  damage  levels  (Figure  5) 
of  the  virgin  material.  (Recall  that  25  different  materials  can  be  Incorporated 
Into  the  finite  element  formulation).  This  model  will  be  referred  to  as  the 
Multiple  Damage  Level  Model.  The  24  damage  levels  were  defined  on  the  basis  of 
equal  multiple  elastic  secant  modulus  reductions.  The  24th  Is  defined  from  the 
ultimate  stress  and  strain,  ou  and  eu,  of  the  true  stress-true  strain  curve. 

Since  each  of  the  24  damage  levels  Is  attained  at  a  different  value  of  strain  en¬ 
ergy  density,  the  true  stress-true  strain  curve  must  be  explicitly  given  (Figure 
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Figure  3  -  Single  damage  level  model  No.  1  uniaxial  tensile  specimen 
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Figure  4  -  Single  damage  level  model  No.  2  uniaxial  tensile  specimen 
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Figure  5  -  Multiple  damage  level  model  uniaxial  tensile  specimen 


5).  For  this  analysis,  the  form  used  was 


a  =  Ac' 


(20a) 


The  constants  were  evaluated  for  a  structural  steel  with  the  following  properties: 


Undamaged  Modulus  E1  =  30.0  x  10°  psi 


Yield  Stress 


Oy  =  35,000  psi 


Ultimate  Stress  ou  =  41,500  psi 

Ultimate  Strain  eu  =  0.0023 


For  this  material,  the  stress-strain  relation  is 


a  *  190500  e 


0.251 


(20b) 


For  the  onset  of  damage  in  the  material,  the  strain  energy  density  in  the  ith 
element  must  be 


<$>,  >  $>ys  <2’) 

For  the  kth  level  of  damage  in  the  material,  characterized  by  an  elastic  modulus 
of  E^,  the  strain  energy  density  in  the  1th  element  must  be 


/dWx  ,dWx  ^  /dW» 

‘  H  i  <w»w 


where,  for  the  form  of  the  given  stress-strain  relation 


A - 


ek  ey  ek 
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The  load  P  Is  the  Independent  variable  In  this  model,  since  all  elements 
have  the  potential  to  sustain  damage  at  a  given  load.  The  Increment  of  the  load 
P  can  be  arbitrarily  chosen  within  the  limits  of  upper  and  lower  bounds  (See  Ap¬ 
pendix  A).  This  Is  In  contrast  to  the  Single  Damage  Level  model,  where  the  dam¬ 
age  level  Is  the  Independent  quantity  governing  the  applied  load  (Figures  6  and 
7). 

The  concept  of  the  mean  damage  level  ¥  and  damage  center  (xd,yd)  have  been  In¬ 
troduced  In  order  to  quantify  the  damage  zone  In  terms  of  a  scalar  magnitude  and 
Its  position  In  the  panel  specimen.  The  mean  damage  level  if  Is  defined  on  the 
basis  of  the  discrete  reductions  In  elastic  modulus.  Defining  the  damage  In  the 
1th  element  as  the  fractional  reduction  in  the  elastic  modules  by 
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The  mean  damage  level  <T  Is  calculated  as 


-  ^  fd1vi 


Elements  which  have  not  reached  the  (dW/dV)y  damage  threshold,  of  course,  have 
values  of  dj  equal  to  zero.  The  coordinates  0Td,yd)  of  the  damage  center  are  de- 


Figure  6  -  Single  damage  level  model  Figure  7  -  Multiple  damage  level  model 


fined  on  the  basis  of  the  first  moment  of  the  damage  distribution 


(26a) 


(26b) 


where  nd  Is  the  number  of  elements  having  damage  at  the  applied  load.  This  Is 
simply  the  discretized  form  of  first  moment  Integrals  with  the  density  function 
taken  to  be  the  damage  level.  The  and  y^  coordinates  are  the  coordinates  of 
the  1th  element  centroid. 


The  locus  of  the  damage  center  coordinates  at  each  load  serves  as  an  Indica¬ 
tion  of  how  the  material  damage,  as  defined  by  the  models  Incorporated  here.  Is 
being  distributed  in  the  center  cracked  panel.  The  crack  tip  Is  the  origin  of 
the  locus  of  points,  and  the  subsequent  curve  Is  Indicative  as  to  whether  the 
lower  levels  of  damage  at  the  periphery  of  the  damage  zone,  or  the  higher  levels 
nearer  to  the  crack  tip  are  predominating  the  damage  accumulation  process  at  a 
given  load  Increment.  Also,  the  effect  of  the  free  boundary  ahead  of  the  crack 
can  be  evaluated  as  to  Its  ability  to  Influence  the  damage  zone's  direction  of 
propagation.  Typically,  one  would  expect  the  damage  zone  boundary  to  propagate 
more  toward  the  free  surface  as  the  zone  size  becomes  significant  with  respect  to 
the  uncracked  ligament  length. 


LOAD-DISPLACEMENT  RESPONSE 


Two  load-displacement  curves  are  given  that  were  generated  from  the  Single 
Damage  Level  model  (Figures  8  and  9).  The  application  of  these  models,  with 
"damaged"  modulus  values  of  0.8  and  0.6  times  the  undamaged  modulus  respectively, 
results  in  decreasing  stiffness  values  of  the  center  cracked  panel  (Figure  11). 
Qualitatively,  the  nonlinear  global  behavior  of  the  panel  specimen  has  been  re¬ 
covered  by  the  application  of  a  simplified  model  accounting  for  material  damage 
in  the  form  of  microcracks. 

One  load-displacement  curve  Is  given  which  was  generated  from  the  Multiple 
Damage  Level  model  (Figure  10).  Recall  that  the  25th  set  of  material  properties 
of  the  Multiple  Damage  Level  model  are  Identical  to  those  of  the  damaged  material 
properties  used  In  the  second  case  of  the  Single  Damage  Level  model.  A  comparison 
of  this  nonlinear  curve  to  those  already  mentioned  shows  that  the  Multiple  Dam¬ 
age  Level  model  results  In  a  load-displacement  curve  and  panel  stiffness  similar 
to  the  first  Single  Damage  Level  model.  This  Is  to  be  expected  since  the  mean 
secant  elastic  modulus  In  the  damage  zone  Is  found  to  be  between  2.5  x  10  and 
2.6  x  10'  psl,  greater  than  either  case  of  the  single  level  model  but  closer  to 
the  first  Single  Damage  Level  model  (Figure  12).  This  Indicates  that  both  models, 
one  based  on  damage  Increments  and  one  on  load  Increments,  are  capable  of  producing 
reasonable  global  response  estimates  of  the  panel  specimen  response  to  applied 
loads. 

The  propagation  of  the  damage  zones  (Figures  13,14,15)  Is  seen  to  be  somewhat 
different  In  each  of  the  cases  considered.  The  shape  of  the  damage  zones  for 
the  Single  Damage  Level  models  shows  significantly  more  damage  directly  ahead  of 
the  crack  for  the  second  case.  This  fact,  which  results  In  a  smaller  undamaged 
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Figure  9  -  Single  damage  level  model  No.  2  load-displacement  curve 
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Figure  11  -  Center  cracked  panel  secant  stiffness  versus  applied  load 
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Figure  12  -  Average  secant  modulus  versus  applied  load 
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ligament  length,  is  probably  the  primary  cause  of  the  decreased  global  stiffness 
of  the  second  case.  Thus,  higher  damage  levels  in  the  enclaves  on  each  side  of 
the  crack  tip  tend  to  enhance  the  strain  energy  density  field  ahead  of  the  crack 
tip.  This  produces  damage  ahead  of  the  crack  at  lower  applied  loads.  For  lower 
levels  of  damage,  this  effect  seems  to  be  reduced,  preserving  the  undamaged  liga¬ 
ment  at  higher  applied  loads. 

Despite  the  differences  In  damage  zone  propagation,  the  damage  center  loci 
(Figure  16)  show  interesting  similarities.  The  shape  of  all  three  loci  are 
basically  the  same.  Initially,  for  loads  less  than  1.8  x  1C5  pounds,  the  damage 
center  follows  the  direction  eQ  =  l-2v  of  maximum  strain  density  in  the  neighbor¬ 
hood  of  the  crack  tip,  as  discussed  by  Sih  [13]  for  a  linear  elastic  body.  At 
higher  values  of  the  applied  load,  the  damage  center  begins  to  move  toward  the 
free  surface  ahead  of  the  crack  tip.  This  seems  to  indicate  that  the  effect 
manifests  itself  at  a  particular  load.  Independent  of  the  specifics  of  the  damage 
models  being  used.  This  movement  of  the  damage  center  may  be  significant  In  re¬ 
lation  to  the  onset  of  macrocrack  growth.  It  may  be  produced  as  a  result  of  the 
damage  zone  size  and  geometry  Itself,  the  Influence  of  the  free  surface  ahead  of 
the  crack  tip,  or  a  combination  of  the  above. 

The  mean  secant  elastic  modulus  (Figure  12)  remains  relatively  constant  with 
respect  to  the  applied  load  for  the  Multiple  Damage  Level  model.  (This  is,  of 
course,  by  definition  the  case  for  the  Single  Damage  Level  model).  Thus,  the 
net  increases  in  damage  at  higher  loads  near  the  crack  tip  are  offset.  In  the 
averaging  process,  by  the  comparatively  lower  damage  levels  developing  In  the 
peripheral  areas  of  the  damage  zone. 


The  suppression  of  macrocrack  growth  as  applied  In  these  analyses  undoubted¬ 
ly  has  an  effect  on  the  damage  zone  development  In  the  panel.  However,  the  con¬ 
sideration  of  this  artificially  constrained  situation  is  essential  In  order  to 
separate  and  explain  the  combined  phenomena  of  micro  and  macro  damage  propagation 
In  the  panel  specimen.  The  use  of  strain  energy  density  In  the  criterion,  as  ap¬ 
plied  here,  describes  the  level  of  damage  In  the  material  around  the  crack  tip. 

It  can  be  applied  with  consistency  to  the  propagation  of  the  macrocrack  Itself, 

In  accordance  with  the  theory  presented  by  Slh  [12,13],  thus  unifying  the  criteria 
for  irreversible  processes,  of  the  types  considered,  at  two  different  scale  lev¬ 
els.  It  should  be  recalled  that  the  microcracks  In  reality  are  accompanied  by 
microplasticity  In  the  sense  of  dislocation  movement  and  the  development  of  slip 
lines  near  the  microcrack  tips.  This  dissipative  mechanism  has  not  been  Included, 
thus  keeping  the  problem  elastic  between  each  damage  or  load  increment,  and  dissi¬ 
pation  only  In  the  sense  of  microcrack  formation  (with  macrocrack  extension 
treated  In  the  second  part  of  this  work).  This  Is  In  contrast  to  the  more  com¬ 
mon  approach  of  using  the  von  Mlses  criteria  of  classical  plasticity  theory  to 
describe  the  size  and  shape  of  the  plastic  zone  around  the  crack,  and  a  separate 
crack  propagation  criteria  based  on  either  a  maximum  stress  component,  strain 
component  or  stress  Intensity  factor. 


CONCLUDING  REMARKS 


The  application  of  two  different  damage  accumulation  models,  based  on  the  . 
strain  energy  density  level  In  the  material  at  the  continuum  scale, -has  been" used 
to  generate  the  global  load-displacement  curves  for  a  center  cracked  panel  In 
the  absence  of  macrocrack  propagation.  These  models  have  qualitatively  shown  v 
how  such  a  criterion  has  the  potential  to  be  Interpreted  In  terms  of  a  specific 
damage  mechanism  at,  or  below,  the  continuum  scale.  Further  Investigation  should 
address  the  details  of  these  damage  mechanisms  end  the  way  lit  which  they  Influ*’* 
ence  the  effective  material  properties  of  the  material.  This,  In  turn,  would 
lead  to  Improved.damage  accumulation  models  end  a.  better  understanding  of  the  . 
effects  material  damage  has  on  the  behaylor  of  structural  members.  - 

Damage  models  which  Incorporate  two  or  more  mechanisms  undoubtedly  would  re¬ 
sult  In  a  more  complete  formulation  of  the  problem.  The  details  of  initiation- 
and  propagation  processes  of  microcracks  will  provide  Insight  Into  the  energy  „ 
dissipating  processes  which  occur  at  this  scale.  Kith  this  Information,  the  ten¬ 
sile  test  specimen,  dimensions  may  be  specifically  chosen  on  the  basis  of  the  dam¬ 
age  processes  which  the  specimen  seeks  to  reflect.  The  continuum  properties  ob¬ 
tained  from  the  specimen's  true  stress-true  strain  curve  will  then  correspond  to 
the  processes  that  are  actually  occurring  In  a  structural  member.  The  appropriate 
choice  of  variables,  which  will  be  capable  of  describing  the  state  of  damage  In 
the  material,  can  then  be  Incorporated  Into  the  constitutive  relations  for  the 
material.  These  relations  will  directly  relate  material  damage  at  the  continuum 
scale  to  the  global  behavior  of  a  structural  member. 
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APPENDIX  A  -  DISCUSSION  OF  UPPER  AND  LOWER  BOUNDS  OF  GENERATED 
LOAD  DISPLACEMENT  CURVES 

In  the  Single  Damage  Level  Model,  the  sequence  In  which  the  finite  elements 
of  the  center  cracked  panel  specimen  fall  Is  unique,  since  the  value  of  the 
strain  energy  density  In  each  element  Is  the  result  of  a  unique  stress  field. 

When  multiple  element  failure  Is  used  to  generate  the  damage  Increment  In  this 
model,  the  question  arises  as  to  how  much  of  a  damage  Increment  can  be  Introduced 
while  still  producing  a  physically  realistic  load-displacement  curve.  Discussion 
of  reasonable  upper  and  lower  bounding  procedures  Is  Included  here  for  this  pur¬ 
pose. 

Consider  the  mth  portion  of  the  generated  load  displacement  curve  (Figure 
A1.1)  for  the  center  cracked  panel  with  a  stiffness  Sffl.  Prior  to  accumulating 
the  damage  resulting  In  this  stiffness,  which  occurred  at  an  applied  load  of 
Pm-l*  the  st1ffness  of  was  S^-,  for  the  m-1  segment  of  the  curve. 

Note  that  in  the  context  of  the  pseudo-linear  elastic  material  formulation  used 
In  the  damage  model,  all  portions  of  the  load-displacement  curve  (except  those 
corresponding  to  the  damage  Increment)  extrapolate  back  to  the  origin  of  the 
P-6  axes.  Regressing  further  along  the  curve  another  Increment  yields  the  m-2 
segment,  with  local  stiffness  and  the  applied  load  P^* 

Define  the  tangent  stiffness  S*  as  the  slope  of  the  line  passing  through  the 
points  on  the  Incremental  curve  labelled  "A”  and  "B".  Thus, 

s* .  Vrfcz  (ai) 

m  6b-6a 
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Figure  Al.l  -  Upper  and  lower  bounds  for  load  and  displacement  Increments 


This  Is,  of  course,  simply  an  approximation  to  the  actual  tangent  stiffness 
3P/36. 


From  the  physical  behavior  of  a  center  cracked  panel  made  of  a  ductile  mate¬ 
rial,  In  the  absence  of  crack  propagation,  one  would  expect  a  smooth  load-dis¬ 
placement  curve  to  develop.  The  slope  of  the  curve  should  remain  constant  or 
decrease.  Implying 


1_  (IE)  <  o 

36  '36'  -  u 


(A2) 


This  condition,  as  applied  to  the  Incremental  load  displacement  curve  through 
the  value  of  S*,  sets  an  upper  bound  to  the  applied  load  Pm. 

Denote  by  a  dark  circle  V  those  (P,6)j  points  on  the  mth  local  stiffness 
curve  where  the  jth  element  will  fall.  These  elements  are  below  the  upper  bound 
set  by  the  Intersection  of  the  S*  and  the  S  lines  which  define  an  upper  bound 
applied  load  pupper*  All  other  elements,  labelled  with  an  open  circle  Mo",  which 
have  not  accumulated  damage  at  an  applied  load  less  than  or  equal  to  P^-j  must 
lie  on  the  secant  modulus  line  at  points  corresponding  to  loads  Pj  which  lie 
above  P^^.  Failure  of  any  of  these  elements  would  produce  a  physically  un¬ 
realistic  load-displacement  curve. 

Alternately,  a  tangent  stiffness  S*  may  be  defined  as  the  slope  of  the  line 
passing  through  the  points  A*  and  B’  on  the  Incremental  load-displacement  curve 
(Figure  A1.2).  In  this  case 


This,  also.  Is  an  approximation  to  the  actual  tangent  stiffness,  dP/dS. 

If  the  physical  behavior  of  the  center  cracked  panel  Is  considered.  It  should 
be  clear  that  due  to  Identical  reasons  as  discussed  above. 


(A4) 


l.e.,  the  tangent  modulus  should  not  Increase  with  Increasing  load.  This  condi¬ 
tion  will  be  satisfied  If  the  nth  damage  Increment  produces  a  displacement  Incre¬ 
ment  A£m  such  that  the  displacement  subsequent  to  the  damage  accumulation  Incre¬ 
ment  Is  greater  than  that  defined  by  the  tangent  stiffness  line.  The  failure  of 
an  undamaged  element,  (along  with  those  undamaged  elements  having  a  higher  strain 
energy  density  for  the  applied  failure  load).  Is  denoted  by  a  dark  circle  On 

the  mth  local  stiffness  curve,  this  defines  the  groups  of  elements  which  satisfy 
this  lower  bound  by  producing  a  sufficiently  large  displacement  Increment  A«m. 

The  other  undamaged  elements,  labelled  with  an  open  circle  "o",  define  those 
sets  of  damaged  elements  which  would  not  satisfy  the  lower  bound  since  an  Insuf¬ 
ficient  displacement  Increment  AAm  would  be  produced. 

In  the  case  of  the  Multiple  Damage  Level  model,  the  same  approach  Is  used 
with  the  following  difference.  The  load  Increment  Is  the  Independent  variable, 
requiring  that  any  load  Pm  satisfying  the  condition  <.  Pffl  <_  satisfies 

the  upper  bound.  The  resulting  damage  Increment  A6m  must  In  turn  satisfy  the 
lower  bound. 

In  application,  it  should  be  noted  that  the  upper  bound  can  be  used  to  Insure 
that  the  applied  load  P^^  Is  not  surpassed  In  choosing  the  damage  Increment. 
The  lower  bound,  however,  can  only  be  used  to  check  on  the  choice  of  the  applied 
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load  P  ,  since  the  displacement  Increment  must  be  calculated  In  order  to  be  com- 
m 

pared  with  the  lower  bound  tangent  stiffness  line. 
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Damage  Increment  Damage  Increment 


APPENDIX  2.1  -  SINGLE  DAMAGE  LEVEL  MODEL  NO.  1  ELEMENT  DAMAGE  SEQUENCE 


No.  1 

56 

57 

55 

66 

No.  2 

67 

73 

44 

No.  3 

45 

62 

52 

No.  4 

68 

74 

36 

37 

No.  5 

43 

46 

63 

38 

Element  Damage  Sequence  (First  -*•  Last) 

55  66  72  58  60  59  53 


WENPIX  2.2  -  SINGLE  DAMAGE  LEVEL  MODEL  NO.  2  ELEMENT  DAMAGE  SEQUENCE 


Element  Damage  Sequence  (First  Last) 

55  66  72  58  60  59  53 


No.  1 

56 

57 

55 

66 

72 

58 

No.  2 

67 

73 

44 

No.  3 

62 

45 

68 

74 

52 

No.  4 

63 

46 

69 

75 

38 

43 

No.  5 

37 

36 

47 

64 

70 

76 

76  39  29 
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APPENDIX  2.3  -  MULTIPLE  DAMAGE  LEVEL  MODEL  ELEMENT  DAMAGE  LEVEL  PRODUCED  BY  APPLIED  LOAD  2 
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